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Abstract 

We address the nonperturbative effects associated with soft charm quarks in inclusive B 
decays. Such corrections are allowed by the OPE, but have largely escaped attention so 
far. The related four-quark 'double heavy' expectation values of the form {B\bccb\B) are 
computed in the 1 / mc expansion by integrating out the charm field to one and two loops. 
A significant enhancement of the two-loop coefficients is noted. A method is suggested 
for evaluating the expectation values of the higher-order 6-quark operators required to 
calculate charm expectation values, free from the overly large ambiguities of dimensional 
analyses. The soft-charm effects were found generally at the level of 0.5% in rsi(6— >-c); 
our literal estimate is somewhat smaller as a result of partial cancellations. We propose 
a direct way to search for such effects in the data. Finally, we discuss the relation of the 
soft charm corrections in inclusive decays with the 'Intrinsic Charm' ansatz. 



* On leave of absence from St. Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188300, Russia 



1 Introduction 



Inclusive weak or radiatively induced electroweak decay rates of heavy flavor, in particular 
of beauty hadrons can be rigorously treated in the heavy quark expansion employing a 
local operator expansion in powers of l/mg (OPE). For fully integrated rates there are 
no nonpcrturbativc corrections scaling like l/mq [1,2], and the leading contributions are 
given to order 1 / ml in terms of the expectation values of the kinetic and chromomagnetic 
heavy quark operators. They incorporate effects due to bound-state dynamics in the initial 
state as well as the propagation of the hard final state quarks through the nonperturbative 
hadronic medium during the decay. In this respect they are analogous to the leading power 
correction in the vacuum current correlator due to the gluon condensate oc [3] . 

At order 1/ml four-quark operators of the form {B\brq qrb\B) appear in the expansion 
of the decay rates. They describe, for instance, weak annihilation (or weak scattering in 
baryons) and Pauli interference which distinguish the decays of hadrons with different 
fiavors of the light spectator quarks - yet they generate fiavor-singlet effects as well. 
Their meaning is transparent: the propagation of soft light quarks in the decay is not 
perturbative, neither literally nor in the duality sense, and the overall contribution from 
these kinematics should be regarded as an additional nonperturbative input [4,5]. It 
parallels the quark condensate effect in the vacuum correlator [3]. 

Inclusive semileptonic B decays offer the cleanest theoretical environment for precision 
application of the OPE. They allow to extract both the heavy quark parameters m^, rric, 
/i^ etc. and \Vcb\ or \Vub\ with (almost) unmatched ^ accuracy (see, for instance, Ref. [7]). 
It is then essential to analyze all contributions, even small ones [8]. 

The four-fermion operator 67^(1— 75) ti'U7'^( 1—75)6 explicitly appears at tree level in the 
semileptonic 6— >it decays [5]. In the dominant semileptonic mode 6— >c analogous contri- 
butions from 67^(1— 75)007*^(1— 75)6 are usually neglected under the tacit assumption that 
c quarks are sufficiently heavy as not to induce nonperturbative effects.^ This assumption 
is valid only up to a point: nonperturbative charm effects could a priori become relevant 
at the one percent level. 

There is a complementary reason to introduce explicitly the four-quark expectation 
values with the charm quark fields, which becomes manifest once the heavy quark expan- 
sion in b^ciu is extended beyond order 1/mf. For instance, in calculating Ftot we expand 
the charm quark Green function in the external field. Similar to differentiating over m^, 
this produces increasingly infrared-sensitive integrals converging at kc ^ rric rather than 
at the hard scale of the energy release. Consequently, the resulting expansion in general 
includes terms scaling like 

To/^hiir— fc— r (1) 



c 



^The b quark mass has been extracted from the threshold cross section of e+e~ bb with similar 
accuracy [6]; that value is in a very good agreement with the recent results from inclusive B decays. 

partial exception is the c-quark-induced Darwin term which, by the equations of motions, reduces 
to the flavor-singlet 47rQs^7oi"^ ^7oi"? operator, manifesting the perturbative mixing of 67^ (1-75)0 07'° (1— 
75)6 with the latter. 
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with / > for k>3. Just for the above reasons these terms describe precisely the lump 
effect of the expectation values like {B\brccrb\B). Including the latter explicitly allows 
to preserve the principal advantage of the OPE for the inclusive rates, viz. expansion in 
l/rub (or, more generally, in the inverse energy release) free from corrections suppressed 
only by the inverse charm mass. 

Therefore, with the charm quark significantly lighter than the b quark and only 
marginally heavy on the QCD mass scale, it is advantageous to explicitly introduce the 
four-quark operators with charm into the OPE and analyze their impact. 

The remainder of the paper will be organized as follows. In Sect. 2 the framework is 
set up which is then used to compute the charm expectation values in the l/rric expansion 
to one- and two- loop level. The expectation values of the 6-quark operators of dimension 7 
and 8 are estimated, and the contributions to rsi(6— >c) from both and CKs/rric terms 
are evaluated. In Sect. 3 we discuss the phenomenology of the soft-charm effects and 
point out that the g^-moments of the decay distributions are rather sensitive to them. 
In Sect. 4 we summarize our conclusions and dwell on a possible connection of these 
effects with 'Intrinsic Charm' in B particles. Appendices guide through technical details 
of our analysis: the conventions are described in Appendix A, the master technique for 
calculation is briefly reviewed in Appendix B and in Appendix C we describe the main 
steps in the technically challenging two-loop calculations. 



2 The framework 



The total semileptonic 6— >c decay rate in the OPE is given by 



r,i(6^c) 



192 7r3 



zo{r){l+Ar\r;^)) 



2m2(/x) 



J 



The coefficients A^^'^^ stand for perturbative ag corrections and r = {mc/mb)" 



ir + 8r — r 



3 ^4 i2r2lnr 



J/ N oi 34 32 „ 2 32 3 10 4 

a(r) — Slnr-\ r — 8r H r r . 

^ ^ 3 3 3 3 



(2) 



(3) 



A detailed discussion can be found in Ref. [8]. We have adopted the Wilsonian prescrip- 
tion of introducing an auxiliary scale /i to separate large- and short-distance dynamics 
consistently. 
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The last term in Eq. (2) proportional to Fg denotes the effect of generic S'?7(3)-singlet 
four-quark operators (other than Darwin operator) of the form hVhqVq with the sum over 
q = u, d, s, and F including both color and Lorentz matrices (to the leading order in 1/mb 
only 7o x 7° or 7^75 x 7*75 structures survive, but one does not need to rely on this). 
Their Wilson coefficients are 0{as), and we neglect these contributions. 

Since the b quark is sufficiently heavy, we consider terms only through order l/mf. This 
includes the four-quark operators with the charm fields, yet only those without derivatives 

- the latter appear to order 1/m^ or higher. Without perturbative loop corrections only 

He = ^{B\hA^-j,)ccY{l-l^)b\B), {-S;+v''v,) , = ^ (4) 

will contribute. This form literally implies the demarcation point /i to be above the c 
quark mass, in which case charm is fully dynamical. The four-quark expectation values 

then are dominated by the large perturbative piece oc fi^/Gn'^. It is natural to eliminate it 
by evolving the theory down to smaller fi and eventually integrating out the charm field. 
This procedure adds nonperturbative corrections to be analyzed below. 

Limiting ourselves to the third order in l/rrib - the leading one for these operators 

- amounts to treating them in the static limit — > 00. There are then in general 
four 0(3)-invariant parity-conserving operators possible of the form bbc'joc and babc^'y^c, 
with two color contraction schemes each. In the actual semileptonic decays a certain 
combination enters which can be decomposed into color-singlet and color-octet operators 
in the s-channel: 

Hc^iOl) = -^{B\H{l-l,)ccY{l-l,)b\B) 

Fc = {Ol) = -^(£|6f7,(l-75)ccfy(l-75)6|5). (5) 

At a high normalization point only He enters at tree level according to Eq. (2), however 
strictly speaking these two mix under renormalization. For our analysis we consider, in 
general, the four operators 



Oy = bb c-)qC 0\ = bab c^^c 

^b c^7oc = b^ab 

which emerge after Fierz reordering of the operators in Eq. (5): 



0°y = Wb c^^oc 0\ = Wab c^ii^c, (6) 



01 - - ^)0t. + 1(1- + ^O^y - ^O", . (7) 

We therefore have vector and axial currents of c quarks created at the origin (the position 
of the b quark); both can be color singlet or octet. 
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For sufficiently lieavy cliarm quarks we can integrate tliem out from the theory and 
have an expansion of Hf. and Ff. (or, in general, for all the four operators in Eq. (7)) in 
terms of the local b quark operators Ok containing only gluon and light quark fields: 



where dk> 3 denote the dimension of the operators Ofe. This series is an expansion in 
inverse powers of charm quark mass. 

2.1 Integrating out the charm quarks 

To integrate out charm quarks from the operators in Eq. (6) we calculate the one-loop 
as well as the two-loop (one hard gluon exchange) contributions. The latter are of the 
leading power in l/nic, however contain an extra power of the running coupling as- 

2.1.1 One- loop expansion 

The operators in Eqs. (6) can generically be written as 



to one loop order the dependence on ri^2 totally factorizes when the charm field is inte- 
grated out. Without hard gluon corrections the expansion (8) reduces to the charm quark 
current in an external gluon field. In the notation of Eq. (9) 



where the average denotes the gluonic operator obtained from the charm Green function 
in the given external gauge field A. The required expectation values are then given by 
the averages of these composite operators over the B meson states. 

In practice we use the well-elaborated Fock-Schwinger gauge method, which yields the 
result directly in terms of gauge- covariant objects. Here we shall only state the results. ^ 
The method is briefly described in Appendix B and discussed in more detail in [10]. 

The absence of the b quark when calculating {cr2c)A brings in additional symmetries 
and simplifies the possible structure of the operators in the expansion. Lorentz and gauge 
invariance ensure that the lowest possible operators are of the general form DGG, with 
D the covariant derivatives and G the gluon field strength. Therefore these operators are 
of order They agree with the results presented in Ref. [9]. The vector color-singlet 

current emerges only at the higher order 

In what follows we will mostly use the matrix representation; the covariant derivative 
Dfj^ = dfj_—iA^ then acts as a commutator. This sometimes will be assumed and not written 
explicitly. More definitions can be found in Appendix A. 

^While the present study was being written a paper [9] appeared which also considered, in a different 
context, the generic charm quark loop in the gluon background to order l/m^. We are grateful to S. Trine 
for bringing it to our attention. Comparing with the expressions there helped us to find an inconsistency 
in the routine used for transforming the results to various forms which lead to incorrect final expressions 
for the 1 /m^ terms presented in the preliminary version of this paper. 





Ji2 = bVib cEsc ; 



(9) 



Ji2 = 6ri(cr2c)A6; 



(10) 



4 



Axial current 

For the axial current with arbitrary color we get to order 



487r^m^ 

where a and b are fundamental (spinor) indices in the SU{3) color space. Note that 
[Da,G°'^] in the first term reduces to the color quark current, glYl/qQlpQ in the adjoint 
representation. For a spatial Lorentz index i we can write 

{A)ba = (2{[A, Bki Ek} + {[A, Eki Bk) - {[A, Bl Ek) + {[A, EklB^}- 

eiu[D^,EkiEi}\^ . (12) 

The application to the semileptonic decays does not require singling out the color- 
singlet piece of the current, and the obtained general expression appears most convenient. 
The case of the color singlet axial current corresponds to taking the trace over the color 
indices; it results in a much simpler equation: 

Al = T,[A] = {c^7>.75C^)a = da Tr G^^G,. (13) 

The axial singlet current has been calculated in Ref. [11]. Assuming that = i is a spatial 
index, the singlet current can be rewritten in the form 

da Tr G'^^G,, = 2diTi (EB) - 4 (Tr EkB, - Tr E^Bk) . ( 14) 

The expression for the general axial current can be cross-checked [10] via the non-Abelian 
axial anomaly using the expansion of the pseudoscalar density. 

Vector current 

For the most general color vector current we get 

-4^[Z}^[D°,[D^,GJ]])^^ . (15) 

It is readily verified that this current is covariantly conserved, D^V^ = Q [10] as required 
by gauge invariance. Note that we have dropped the leading term D^G^fi = gl J2q 

with a coefficient ocln since it is accounted for in the Darwin operator in the width 

at order The pure octet component of the current can be obtained by varying the 

b c 

gauge field effective action with respect to which yields Tr(V^t''); this is explicitly 
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verified in [10]. By setting /x = - the only component surviving the static hmit of the b 
quarks - we can write the expression for the current as 

{Vo)ba = 13(D-(ExB) + D-{BxE))+M(2D,D'' - D,DP)DpEk) (16) 

(here the usual commutators are assumed when we symbolically write products of covari- 
ant derivatives). 

The vector color-singlet current vanishes at order as follows from an explicit 

calculation. This was anticipated in Ref. [11] based on C-parity arguments. The leading 
term appears at order 1/m^ and has the general form DGGG. This current describes 
the effective local interaction of the photon with the gluon fields (in the vacuum) due to 
heavy quarks, and is of independent interest. It can be written as a total derivative: 

^ {ei,e)A = - Yii^^ {u{G^pG^^G,,U + HG.pG^'G^,},) , (17) 

where {. . means that the symmetric matrix product ~ dabc of the three matrices 

must be taken. This expression vanishes for the SU{2) gauge group, as expected from 
G-parity [10,11]. Gauge invariance ensures this current to be conserved, d'^V^ = 0, which 
is evident from the symmetry properties of the field strength tensors. The two terms, 
however, do not vanish for the Abelian U{1) group; they can be obtained from the Euler- 
Heisenberg effective action as a variation with respect to the gauge field [10]. Eq. (17) 
therefore represents the minimal non- Abelian extension of the Abelian current. Once 
again anticipating that = in the expectation values, the two terms above can be 
written as 

Vo^ = a"(14{G,^G^^G^o}d + HGj^G^'G^oh)^ 

d,{A{EaE^) - 24(£;„S2) - 7{B,{E ■ B)}^ . (18) 
The one-loop contributions considered above add up to 

Hl-'°°^^~{bVob) + \{baAb). (19) 
2.2 Two-loop expansion to order cx.s{mc) /rric 

Accounting for color Coulomb effects of the b quark reduces the symmetries of the c quark 
interaction with the soft gluons and allows more operators to appear. In particular, gluon 
operators of dimension 4 suppressed only as l/rric can appear for both vector and axial 
currents at the price of an extra loop and the associated factor as{mc)- There are six 
such effective operators allowed on general grounds by P and T parities:^ 

Es = ^S^bbE--E'b Bs = ^SabbB^-B^b 

Ed = IdabcbE'' ■ E%b Bd = \dabcbB'' ■ B%b (20) 

Ej = -iUM" X E'^-atcb Bf = -\UJbB'' x B^'-atJ) . 
^The operators with the chromomagnetic field were left out in Ref. [8]. 
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The first four operators appear in the vector current while the last two are generated 
by the axial current; the former lead to the spin-singlet operator and the latter to the 
operator proportional to the b quark spin. 

Conventional wisdom would tell us that the extra loop and the associated suppression 
by the perturbative factor ^a^/Tr would yield a very small factor. It was conjectured in 
Ref. [81 that the enhancement factor in the matrix element would fall short of com- 

Mhadr 

pensating for it. Nevertheless, we have calculated these contributions and actually found 
a significant enhancement of the perturbative corrections over the naive loop counting 
rules. It is apparently associated with the physics of the static Coulomb field: the typical 
coefficient contains tt^ on top of the usual loop factor ag/n. Therefore these higher order 
contributions may actually dominate unless a cancellation among different operators at 
this order takes place. 

Explicit calculations yield to order as/rric- 

QJ 

{Ov) = 7^ (4708(B,)-1511(£;d) + 9752(^5) - 2470(£;5)) 



21133 TT me 

a 

26325 Trm 

Ols 

2933 Trm, 

The average denotes {X) = ^{B\X\B). Using Eqs. (7) we obtain 
^0''"°" = 7^1^ [74100(^5) + 39810(^rf) - 17720(E/) - 292560(5^) - 144600(5^) 



(Oa) = .,025!^ (-181(g;) - 185(i^;)) 

(01) = ;^^^^(-422(S,)- 295(£;,)) (21) 



- 22672(fi/)] 

-[ -121928(2(^5 + Bd) - Bf) - 189944(5^ + Bd + Bf) 



212345 nrric 

+ 57530{2{Es + E^) - Ef) + 4370(£;5 + Ea + Ef) 

- 273375(^5 - Es) - 170505(^5 + Es) ] . (22) 

The representation above is adapted to the estimates of the matrix elements given in the 
next section. 

2.3 Estimates of the expectation values 

Our goal is to provide realistic estimates of the scale of the expectation values of the opera- 
tors obtained by integrating out the charm field and to get a better idea of the importance 
of the nonperturbative charm effects when they are treated in the l/rric expansion. Esti- 
mating the higher-dimensional expectation values is a notoriously difficult problem, and 
our primary task is to elaborate a method free from arbitrarily appearing huge factors like 
powers of Att which all too often plague meaningfulness of 'Naive Dimensional Analyses'. 
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Since the two-loop induced D = 7 operators have large Wilson coefficients, we start out 
scrutinizing these contributions. We have elaborated more certain estimates of five out 
of six D = 7 operators in Eqs. (20). Our estimates for the one-loop D = 8 and higher 
operators may be less accurate, yet indicate that their effect is presumably small enough 
not to introduce significant uncertainties in rsi(-B). 

2.3.1 Two-loop l/rric expectation values 

The expectation values oi D — 7 operators in Eqs. (20) have been denoted {Eg), {Es), 
respectively (the B* expectation values are related to them in an obvious way, but 
we do not need them). Let us first consider the 'color-through' operators which do not 
have the color trace for gauge fields separately, and we start with those containing the 
chromomagnetic field. In what follows we often pass to the quantum mechanical notations. 
For instance, we assume 

b'Kkni...'Kjb — )• 7rfc7r;...7rj|QM nk ^ -iDk, ttq = iDo-nib . (23) 

Indeed, as a general rule / d^x d^yb7rkb{0,x) b7rib{0,y) — J d^yb7rk7rib{0,y) holds for static 
quarks, since at equal time 6'(0, x)b'{0, y) — S^{x — y)5^-'' {i, j arc color indices), and for the 
quantum mechanical local heavy quark operator O the notation in second quantization 
reads as 

Oqm = / d^xbOb{0,x) . (24) 

The operator ttq acting on the right corresponds to —Ti, or more accurately, it plays the 
role of the commutator of —7i with what follows it. 

We estimate the S-field expectation values using a generalization of the factorization 
ansatz. More precisely, it takes the form of the saturation by the (multiplet of) the ground 
heavy quark states. For a general analysis, it is convenient to consider the static limit, 
where the spin of the heavy quarks gets decoupled. In this world the usual B and B* 
mesons belong to a single spin-| hadron state |0o), see Ref. [12], with 

^2 2 

{no\bnj7ikb\no) = ^Sjk'^l^o - ^^'JcTj-fc^o, <yjk = iejki oi ■ (25) 

Setting 

(nol^i^ikl^^o) - E(^o|5,|^^o)(^ol^fel^o) (26) 

A 

where we have shown explicitly the polarization of the intermediate state, we get 

(no|5,-5,|no) ~ ^^^Ja.a.^o • (27) 
Translating this relation to the world of actual spinor b quarks is simple. For instance, 

— * — * — * 

a — —2j, and in B meson j — —sq, where sq and j are the spin and the angular 
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momentum of the the b quark and of the hght cloud, respectively. With — Sq — ^ and 

— * n 

juQ — — I we have 

{Bs) + {Ba) = {B\bB%\B) ~ ~ 0.041 GeV^ 

{Bf) = {B\bBxB -ablB) ~ ^Ml! ~ 0.082 GeV^ (28) 

This is clearly consistent with the direct factorization for the invariant combination 

{Bs) + (Ba) + {Bf) = {B\b{aBfb\B) ~ {B\baBb\Bf = {filf . (29) 

The combination {2{Bs+Bd)-Bf) = {axB + 2iB){axB -2iB) b\B) vanishes in the 

factorization approximation since {B x j — iB)\Qo) — (in general, j does not commute 
with B, yet its commutator vanishes if projected onto the ground state Qq). Since the two 
linear combinations are squares of the above operators, the factorization estimate gives a 
lower bound for both expectation values. 

A similar ground-state factorization cannot be used for the operators with the chro- 
moelectric field: in the saturation of their product 

{Qo\EjEk\Qo) = ^{no\Ej\n) {n\Ek\no) (30) 

n 

only the P-wave states \n) with parity opposite to that of flo contribute. However, since 
Ek = i[7^o, '^k] we have 

{no\Ek\Pn) = i£„(Ookik|P„), £„ = M^-Ma,. (31) 

Moreover, P„ can be either j = | or j = | states. The ^-states are produced by the 

— * — * — * „ 

(pseudo)scalar product aE; the operator ax E — 2iE creates ^-states. Paralleling the 
derivation of Ref. [12] we use 

{x^"'^\n,\no)=e^^X^^^O: (32) 

where the spinor 0^") and the Rarita-Schwinger spinor x^™-* (obeying aiXi""^ — 0) 
scribe wavefunctions of the j — ^ and J = | excited P-wave states, respectively. For the 
chromoelectric field we then simply get extra factors of 

{Qo\bEjEkb\Qo) = 



£ T, 



(m) 1 2 



3/2 



+ 



E 



e„\T- 



(")|2 
1/2 I 



E 



T, 



»|2 



3/2 



-E ^^1^' 



W|2 
1/2 I 



(33) 
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compared to the sum rules for the products in Eq. (25): 

^ = el\rij^\- + Y: 4|r{/1r, f = 2^ eUrijJl' - ^l^i/ll' ■ (34) 

^ m n ^ rn n 

The result can be anticipated knowing the above decomposition into the |- and |- 
projectors leading to a relation similar to the case of chromomagnetic field operators: 

2{Es+Ed)-Ef ^h{axE+2iE){axE-2iE)b, Eg+Ea+Ef ^b{aE)'^ b . (35) 

The former operator is saturated by |-states and the latter by |-states. The two sum 
rules stemming from Eq. (33) can be written as 

{2{Es+Ea) - Ef) = el/,{2^l + 4) (36) 
{Es + E, + Ef) = el/M-fil) (37) 

where £1/2,3/2 denote the average excitation energies for the two families. 

In a direct analogy to the ground-state factorization ansatz for the chromomagnetic 
field operators given above one retains only the contribution of the lowest P-wave states for 
both |- and |-families, m = n = l in the sum over the intermediates states. This amounts 
to equating the mass gaps e in Eqs. (36), (37) to the lowest corresponding excitation 

energy, £3/2 = £3^2 = £^3/2 and £1/2 = £^/2 = ^1/2 ) and yields a lower bound which is 
simultaneously the factorization estimate. 

The hierarchy for the operators with the chromoelectric field appears to be opposite 
to what we get in case of the chromomagnetic field. The experimental evidence that 
III — /Iq ^ - i.e., the proximity to the so-called 'BPS limit' [13] - implies that the 
combination 2{Es+Ed+Ef) must be particularly suppressed: for in the BPS limit all T1/2 
vanish. The 'BPS' suppression of (Es+Ea+Ef) parallels the smallness of nonfactorizable 
average {2{Bs+Bd)—B f) , while the saturation of {2{Es+Ed) — Ef) by the lowest | P-wave 
excitation is a counterpart of the ground-state factorization for {Bs + B^+Bf). 

It is interesting to note that the saturation of the small velocity (SV, Shifman-Voloshin) 
sum rules by the lowest intermediate state holds with amazing accuracy [14] in the exactly 
solvable 't Hooft model [15]. Without spin in two-dimensional QCD there exists only a 
single family as analogues of the P-wave states. It is intriguing that the latest experiments 
seem to indicate a similar good saturation in actual QCD (for a recent discussion see 
Ref. [16]), although the data so far allow to address this question with a reasonable 
accuracy only for the |-channel. 

The ground-state saturation for the magnetic field hke used in Eq. (26) cannot be 
studied in two-dimensional models in view of the single space dimension. For the sym- 
metric products of the momentum operators, viz. the kinetic operators, the ground-state 
saturation is less accmatc, however the deviation again is almost completely generated 
by the first radial excitation. 

Due to the proximity to the 'BPS' regime, /x^— /x|r<^//^, the estimate of the suppressed 
combination Es+E^+Ef has significant relative uncertainty, yet the absolute uncertainty 
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must be small. The smallness of its coefficient in H^, Eq. (22) makes the uncertainty 
unimportant in practice. For numerical estimates here and in what follows we use the 
following values of the heavy quark parameters A ~ 0.63 GeV, fi^ ~ 0.4 GeV^, ji^ ~ 
0.35 GeV^ pI, ~ 0.2 GeV^ pis ~ -0.15 GeV^ while £3/2 ~ 0.4 GeV and £1/2 ~ 0.35 GeV 
yielding literally 

(2(E5 + Erf)-E/)~0.18GeV^ 

Es + Ed + EfC^ 0.006 GeV^ (38) 

Equating £3/2 and £1/2 results in even simpler relations: 

Es + E^c^ e^l^l, Ef ~ -£Vg ■ (39) 

The operators where the gauge fields E"^ or B'^ are in a color singlet are less certain. 
This is a counterpart of the observed pattern discussed recently in Ref. [17]: factorization 
often fails where vacuum scalar contributions are possible. We nevertheless can estimate 
the combination Bs — Es since it represents (in the chiral limit for light quarks) the 
expectation value of the density at the origin of the trace of the energy-momentum tensor 
O^n of the light degrees of freedom, nonvanishing due to the scale anomaly: 

On the other hand, the space integral of this density over the B meson volume amounts 
to K = MB-mb [18]: 

A=f^(^|TrGj.(0)|5). (41) 
Therefore, by dimensional estimates we would have 

where ^tiRq is the characteristic volume of the B meson bound state. Rq ^ l//ihadr- To 
make the estimates of the effective volume less vulnerable against arbitrariness in the 
powers of 27r inherent in translating energy scales into inverse characteristic distances, we 
use the refined approach based on the useful relation derived in Ref. [19]. Namely, for a 
color singlet hght-field operator A a general relation holds 

{B\hAhmB) = / (0-^a(?) , (43) 

where denotes the transition formfactor of the operator A alone: 

{B{q)\AmB^=,) = :FA{q) . (44) 
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We apply this relation to the energy-momentum trace operator, 




Eq. (41) then fixes J^g^{0) = A. Assuming Ta^iq^) = J^G2(0)e-«'/^' with M^-O.SGeV^ 
and A ~ 0.63 GeV we get 




{M^/^A ~ 0.005 GeV^ (45) 



and 



(Bs-Es) 



IGtt^ 1 



{B\b 



Tr GlXQ)\B) ~ -^V^{M^)^/^A ~ -0.03 GeV^ 



27 2Mb 



87rQ;2 



(46) 



Let us note that the 'vacuum factorization' estimate 



(5166(0)15) (^TrCyo = TV G^o 



is inappropriate since it represents the disconnected piece of the matrix element. For 
this reason it has, e.g., a higher scaling in the number of colors, A^^, while the effects we 
discuss scale only like A^^ ■ We also note that our estimate is negative, necessarily of the 
opposite sign to the vacuum expectation value of the gluon condensate. We interpret this 
as the fact that nonperturbative configurations lowering the vacuum energy are suppressed 
inside B mesons leading to a higher energy density, well in agreement with intuition and 
in the spirit of the assumptions of the traditional bag models for hadrons [20]. On the 
other hand, our estimate (45) taken at face value suggests that the decrease amounts to 
a significant fraction of the vacuum density, while in the framework of the conventional 
1/Nc it would be suppressed as long as the number of colors is large enough. Once again, 
this difference is in line with the picture developed recently and discussed in Ref. [17]. 

We were not able to come up with an equally justified estimate of the B meson 
expectation value of the complementary color-separated contribution bTr{B'^+E'^))b and 



have to rely on general scale estimates of its possible magnitude. We choose to parametrize 
(Bg+Es) in terms of {Bs — Eg), 



Since we speak here about nonperturbative effects and the operator in question is spin- 
independent, we may expect quantitatively that the contribution of the chromoelectric 
field (negative in the vacuum) is suppressed in the B meson, thus increasing (i?^+£'^)(0). 
The quahtative arguments in the spirit of the bag models suggest that (Bs+Es) < K^^ — 
Es) I since they assume the nonperturbative fields to be suppressed in hadrons, therefore 



(Bs+Es) 



h+{Bs — Es) ; 



12 



yielding < < 1. However, there is a significant chromomagnetic field in B meson 
associated with the spin of the light degrees of freedom, ^ 0.35 GeV^, therefore it 
should be admitted that such an argument may not fully apply to the chromomagnetic 
field operators. Assuming naive positivity of the colorless operators Tr and Tr one 
would be led to conclude that h+>l; however, in the case of hadrons with a heavy quark 
the renormalization properties of these operators change, and the naive positivity may be 
lost. 

The tables below summarize our estimates of the D = 7 expectation values. Wc show 
the separate matrix elements for the combinations naturally appearing in our approach 
as described above, and what follows from this for the four operators in Eq. (21). The 
latter numbers assume as{mc) = 0.3 and = 1.15 GeV (a short-distance Euclidean mass is 
adequate here). In Table 2 we also show the estimated contribution to the total B^XcC-v 
decay rate. Taking the numbers at face value we find a strong cancellation among these 
contributions to Fgi, —0.0007(1 + /i+); we consider such a degree of suppression accidental 
and not representative. 



3277^(0^.) 


327r2(0^) 






0.043 - 0.012/1+ 


0.004 + 0.028/1+ 


-0.038 


-0.033 



Table 1: Estimates of the expectation values of the four charm quark current densities at origin, 
at order in GeV^. We assume a5(mc) = 0.3 and mc = 1.15 GeV. 





{Bs+Bd+Bf) 


{Es+Ed+Ef) 


{2{E5+Ed)-Ef) 


{Bs-Es) 


{Bs+Es) 


GeV^ 


0.12 


0.007 


0.18 


-0.03 


0.03/i+ 


<5rsi/rsi 


-0.0034 


5-10-*^ 


0.0015 


0.0012 


-0.00074/1+ 



Table 2: Estimates of the expectation values of the density at origin of the D = A gluon field 
operators, in GeV^, and the corresponding contributions to Ts\{h^c). The expectation value 
of 2{Bs+Bd,) — Bf vanishes in the factorization approximation. 



2.3.2 D — 8 operators at one loop 

Next we estimate contributions from operators generated by integrating out charm at one 
loop which does not incorporate hard gluons with momenta ~mc. As discussed previously 
they are of order I /ml and can be expected at a fraction of percent level. 

Not much is known about the expectation values of higher-dimensional operators. 
Usually one applies dimensional estimates with characteristic hadronic scales entering the 
problems. Although still leaving elements of freedom, they allow to estimate the potential 
magnitude of the effects. The rules for our dimensional estimates are formulated in the 
second part of this section. Unfortunately it is not possible to infer the sign of the 
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contributions in this way. Related to this, we would have to add various contributions - 
which proliferate for higher-dimensional operators - incoherently assuming they have same 
sign. While we consider these counting rules a useful tool for estimating the magnitude 
of various terms in general, they can quite possibly overestimate the net effect on the 
semileptonic width in question. 

On the other hand, we suggest a more credible way to estimate the relevant expectation 
values based on the ground-state factorization. The advantage is that it predicts the sign 
of all the contributions that do not vanish in this approximation, and is free from arbitrary 
factors similar to the conventional vacuum factorization which has been used for decades 
with reasonable success. It appears possible to estimate all the contributions in this way 
since at one loop the effective operators appear directly in the 'color-through' form, cf. 
Eqs. (11), (15), of the generic structure 

hnkT^^-K^-KpTiib , (47) 

with or without the 6-quark spin matrix Bi,. The axial current Aj includes 0"^, while the 
vector current V is spin-singlet. Besides the terms with four spacclike and one timelike 
derivatives coming from the chromoelectric field, both currents contain a term with three 
timelike derivatives. 

For each of the matrix elements with a single ttq there is a unique way to insert the 
ground state \VLq) in the product. When ttq is in the middle, z/ = in Eq. (47), the 
factorization contribution vanishes due to the equation of motion ttq^ = for the h field. 
All the other expectation values can be evaluated using 

Pz) P'ls A^G 

-{QQ\7lk7loni\flQ) = -^<^fe/ + -^CTkl, (QokfcTTilfio) = ~^^kl ^CT^ . (48) 

The estimates of the expectation values of the operators with three time derivatives follows 
from the SV sum rules paralleling the consideration of the D = 7 operators, Eqs. (36-37): 



-(f^okfe7'"o7ro7ro7r/|r2o) = £3/2 — - + £ 



3/2 ^ + ^1/2 ^ I (>ki 



I ~^3/2 ^ ^^1/2 ^ I (^kh (49) 

where £3^2 1/2 average squares of the P-wave excitation energies in the |- and \- 
channels, respectively. One generally has £3^2 1/2 ^^3/2 1/2 compared to the similar average 
energies introduced in the previous subsection, where the lower bounds follow from the 
Holder inequality. Collecting everything and using the same numerical values as given 
above Eq. (38) and £3/2,1/2 — (0.4 GeV)^, we get 

327r2(Vo) ~ -^{l00/.^pl,-35/x^pi5-16(£l/2(2p^-py +£l/2(pl,+py)} 

-0.3GeV^ (50) 
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327r2(a4 ^ " ^ - Gfilpls + l^lpls - 2[^1/2(2p!)-pis) - '^el/^{pl+ pls)]} 

fti -0.2GeV3 (51) 

According to our estimates the vector current dominates 327r^i?p~'°°PRi 0.3 GeV^, Eq. (19); 
the correction to the semileptonic width is then 

^ fti 0.003. (52) 
Tsi i/ml 

As mentioned above, we also tried educated dimensional arguments to estimate the 
higher-order operators. This was done adopting the following rules: 

• The covariant derivative, either spacelike or timelike, counts as 0.35 GeV, and sum- 
mation over the Lorentz index brings in the number of components. For instance, the 
estimate for jji^ would be 0.37 GeV^. 

• Consequently, each factor oi E or B field counts as (0.35 GeV)^ per component. 

• For the anticommutators of the two field strengths we put an extra factor y/2 while 
the commutators of field strengths or covariant derivatives count simply as the product. 

• The spin-dependent operators come contracted with the a matrices. We treat c? as a 
vector with unit components, which according to the above prescriptions reproduces the 
actual chromomagnetic expectation value. 

• All terms in an expression obtained by integrating out charm in a particular current 
are added in moduli regardless of the sign of the corresponding coefficient. 

Using various identities and the equation of motion for the b field, the results for 
higher-order operators can often be cast in different forms resulting in different estimates, 
in particular due to the last rule or due to using the equation of motion for the b field. In 
this case of alternative values we adopt the smallest of the obtained numbers. 

While admittedly of fimited accuracy, these rules seem reasonable, yielding the right 
size for the lowest-dimension expectation values. Applying these to the estimates of D = 8 
operators in the one- loop vector or axial current, we obtain 

327r^(Vo) ~ 0.15 GeV^ 327r^(a^) ~ 0.2 GeV^ (53) 

qualitatively consistent with the factorization estimates; the similarity for the axial cur- 
rent is probably somewhat accidental. The smaller value for the vector current probably 
indicates that counting the commutator with a unit factor in the above rules may under- 
value the higher- dimension local operators. Of course, the dimensional estimates do not 
allow to specify the relative sign of the vector and the axial contributions, and therefore 
would allow a larger overall effect. Adding the two contributions in modulus yields 



^ 3^^2 K rn ? ^^ ^ Q_QQ4 _ ^^^^ 



In the case of the D — 7 operators obtained at two loops we observed significant cancel- 
lations between different terms in the expressions for a particular current. This suggests 
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that the rule of incoherently adding all contributions may noticeably overestimates the 
magnitude of the actual expectation values. That is what we found in estimating the 
matrix elements through the ground-state factorization, although here the difference is 
not large numerically due to dominance of the axial current contribution. 

2.4 Exponential effects 

For a sufficiently heavy quark, in our case charm, the leading operators in the power 
expansion when integrating it out must accurately describe the nonperturbative effects 
associated with these fields. This expansion - like all such power series expansions - is 
as a rule only asymptotic. Therefore it would leave out generic exponential terms scaling 
like e~^™''=/^hadr- ^j-^e exponent, in general, could even be a non-integer power of rric. For 
a dedicated discussion we refer the reader to the papers [21,22]. A related consideration 
actually suggests that twice the charm quark mass enters the exponent. 

It is possible to develop an approach similar to that of Ref. [21] to analyze the asymp- 
totic behavior of the exponential contributions or to estimate their effect in a concrete 
dynamical ansatz for the nonperturbative light fields. We did not attempt this here. 
Nevertheless, charm quarks can be considered as marginally heavy in actual QCD. Vir- 
tual charm effects must definitely be suppressed, but not necessarily to the extent where 
the asymptotic regime can be trusted. They are expected to decrease with growing rric- 
Hence, a safe upper bound would come from similar effects of the nonvalence d quark in 
the semileptonic b^uiu decays of B^* mesons. The latter can potentially be at the scale 
of a couple percent [23]. Assuming the charm mass suppression to be a factor of five or 
stronger, and accounting for the phase space difference, we arrive at the scale 

He ~ 0.001 GeV^ , (55) 

or 5rsi/rsi~0.4%. Whether the truncated power expansion yields such terms, or results 
in smaller effects due to specific cancellations, we do not consider it justified to a priori 
rule out such effects up to the 0.5% level. 

Therefore, the most direct way is to introduce the soft-charm effects as required by 
the OPE, yet described by expectation values taken as free parameters and analyze the 
observable effects they induce. Then one can obtain direct experimental bounds from 
precision measurements of the inclusive observables in 5 — > Xciv decays, thereby placing 
more reliable upper bounds on the limitations they impose on extraction of Vcb- 

3 Phenomenology 

All the inclusive semileptonic B decay distributions are described in terms of five struc- 
ture functions Wi{q'^,qo) representing different Lorentz combinations [24]. In the l/rrib 
expansion through the OPE the effects of the four-quark charm operators appear in the 



16 



structure functions as a delta-function at maximal (f' — {m}) — mc)'^ and qo — mb — mc, 



Swi oc S'^{qf^-{mb-mc)Vi^) oc (5(g^-(m6-mc)^) , S{qo-{mb-mc)) . 



ql-q^ 

The strong interactions smear out this distribution over the range ~ /ihadr ~ Aqcd in 
and in go- However in the fully inclusive characteristics this is an effect of a higher 
order in 1/mb, and we neglect it here. Likewise, we disregard the perturbatively-induced 
tails towards lower values of q^ and qo, associated with the anomalous dimension of the 
four-quark operators [5]. 

Of the five weak decay structure functions of B meson three contribute to the decays 
with massless leptons, 



d^r 



e[q,-E,-^)e{E,)e{q^) X 

{2gV + [AEi,{qQ-Et)-q^]w2 + 2q^{2Ei-qQ)w^} . (56) 



Yet only two of five Lorentz structures, —(5^,/ and v^jVi, are independent at the point 
where q^ — q^ (i.e., g^— >0). These are described by the structure functions Wi and ^2, 
respectively. Contributions from either vanish being proportional to q or reduce to 
that of Wi,2- Two kinematic step-functions in Eq. (56) effectively yield a 5{E(^ 



mt,-mc \ 
2 ) 

type lepton spectrum, as expected (it is also smeared in reality over an interval ~ A*hadr) • 
Furthermore, for massless leptons the effect of W2 vanishes as well, since one has 

rsi = ^f^/^'"'dg2£j^dgoy^(?V(?^,?o) + ^(%'-0^^ ■ 
Therefore, for the correction to the total width in our case we have 

G%\Vrh\^ 



with 



5xrsi = / / ^'\f^^xMQ' Qo) , (58) 

S^Wi{q'^,qo)^^HcS{q'^-{mb-mcf) ^ S{qo-{mb-mc)) . (59) 

Eqs. (56), (58), (59) allow one to evaluate the soft-charm effects in various inclusive 
moments used to experimentally determine the heavy quark parameters, in terms of H^- 
The dimensionless ratio 

X=^4^^ 32.^^^^^^^ 3GeV- • 

determining the correction to the total semileptonic width, sets the overall scale of the 
effects. According to our estimates, it can hardly exceed the one percent level, which 
places the first benchmark for the accuracy required to detect or constrain them further. 
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Since X<^1, in what follows we will keep only the terms linear in x and drop everything 
that scales like higher. 

The average charged lepton energy receives a correction proportional to 

5^{E,)^R^'x{^-{E,)o) , (60) 

where (-E'^)o — 1.38 GeV is the conventional average without the soft-charm effects; it is 
strongly dominated by the parton-level piece. {Rq is the decay fraction left by applying 
kinematic cuts; it!o = 1 for the full moment.) Since the related excess or depletion of the 
rates emerges near ~ Mb-m^* ^ 1.64 GeV which is close to the bare lepton average 
{Ei)o, the impact of the charm field expectation values on (Ei) is strongly suppressed. In 
principle, {Ef )o grows when a lower cut on Ei is applied, while the soft-charm contribution 
remains unaltered up to E^^^^ ~ 1.3 GeV. However, the theoretical precision is probably 
insufficient to detect the mismatch in the overall normalization of the low-E^ part of the 
spectrum at the sub-percent level. 

A similar problem would plague the utility of the higher lepton energy moments: 

^iio^X ((^ - {Ei)o? - {[Ee-{Ei)]\) 
SAEi-{Et)]')^R^'x {{^ - {Ei)or - {[Ei-{Ei)]\) ; (61) 

the numerical values of the moments with the subscript zero are measured in a number 
of experiments and/or are calculated theoretically, see, e.g. Rcf. [25]. The underlying 
feature common to all lepton moments is that they arc dominated by the trcc-lcvcl or 
'partonic' contribution, with all nonperturbative effects being small corrections. At the 
same time the soft-charm effect is not enhanced, but rather suppressed by the smallness 
of the difference between E-^ and (Eq). 

At first glance, the moments of the hadronic invariant mass squared may look better 
since most of the deviation of (Mx) from A'/|), is due to nonperturbative effects, so that the 
'parton background' is suppressed. However, the contributions we focus on are expected 
to populate just the domain of Mx near Mo or slightly above, where the bulk of the 
decays happen: 

s^{m],)c^r^'x{mI,-{m],)o) 

5AMl-{Ml)f)..R^,^x [{Ml. - {MDof - {[Ml-{Ml)]\) 
5AMl-{Ml)f)^R:^'x {{Ml. - {Ml),f - {[Ml-{Ml)]\) . (62) 

Detecting such shifts would require controlling theory predictions for the moments, roughly 
speaking, at a percent level, which does not look realistic, see, e.g. Ref. [25]. (For (M|-) 
this would apply to (M|.)o-M|,* ~0.6 GeV^) 

A realistic way to look for the nonperturbative soft-charm effects is studying directly 
the ^^-distribution, or the associated moments. By their nature the contributions 
in question are located around the maximal (f ~ {Mb — Mjy*)^ ~ 10.7 GeV^, while the 
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average is less than half of this maximal value. Therefore, the moments appear to 
be sensitive: 

5^{q^)c^R^'x{{MB-Mn^f-{q^),) 

<^x([?'-(?')]') ^ ^0 ' X {[{MB-Mr,.f - {q\f - {[q'-{q')]\) , (63) 

where (g^)^ Ri4.3 GeV^ {[q^-{q^)f)Q^l.hGeV^ and {[q^-{q^)f)Q^l GeV^ (without lower 
cut on E(). Similarly, one can study the moments of go instead of the moments of (f. 

In practice, precision measurements of inclusive B decays usually require a lower cut 
on the charge lepton energy around 1 GeV. To facilitate comparison of experimental data 
with theory, we present here the numerical predictions for the q^ moments, which are 
available applying the OPE [25-27].^ For instance, based on the recent fit [7] to the 
existing data on B decays, we get for the lepton energy cut at 1 GeV 

(g2)~ 4.83 GeV2 + 7.2 GeV^X, 
- {<f)f) - 7-53 GeV^ + 33 GeV^x , 
{['f - {(l^)f) - 4-6 GeV^ + 245 GeV^x • (64) 

The dependence on the usual heavy quark parameters can be approximated in the way 
similar to the one adopted in Ref. [25]: 

M{mb, rric, nl, fi^, p%, pls;as) = V + B {mb-A.6 GeV) + C (mc- 1.15 GeV) 

+ P ifil-OA GeV^) + D {p%-0.2 GeV^) (65) 
+ G (/i^-0.35 GeV^) + L {pls+0.15 GeV^) . 

The reference values V and the hnear dependence coefficients B to L are given in Tables 
3 to 5. V has dimension of the moment M. itself and the quoted numbers are in GeV to 
the corresponding power; the values of the coefficients B to L are likewise in the proper 
power of GeV, and Ecut are shown in GeV. We note that the preliminary data reported 
by CLEG are in agreement with theory predictions without significant nonperturbative 
charm effects. Moreover, an estimated accuracy of ±0.1 GeV^ in (g^) and ±0.5 GeV^ in 
translates into the sensitivity in x of 0.015. At the percent level of precision 
the reliability of theory predictions becomes important. Matching the demand would 
require calculating the ag-corrections to the Wilson coefficients of the power-suppressed 

^We give here the resuhs obtained with the same simphfications as in Ref. [25] which proved to be a 
good approximation. The complete cut-dependence in the perturbative corrections is available [26,27]; 
therefore, the estimates can be easily refined in this respect, if necessary. 
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operators, which is presently the hmiting factor. 



Ecut 


V 


B 


C 


P 


D 


G 


L 





4.31992 


2.6176 


-2.1118 


-0.008429 


-1.138 


-0.627 


0.1361 


0.5 


4.41171 


2.6084 


-2.1106 


-0.008009 


-1.147 


-0.6301 


0.1369 


1.0 


4.78331 


2.6724 


-2.2266 


-0.02488 


-1.311 


-0.709 


0.146 


1.5 


5.13019 


3.3376 


-2.9504 


-0.1478 


-2. 


-1.135 


0.1559 



Table 3. First moment of the lepton pair invariant mass {q^). 



Ecut 


V 


B 


C 


P 


D 


G 


L 





7.23244 


9.3514 


-7.2682 


-0.04725 


-7.737 


-3.167 


0.6866 


0.5 


7.1329 


9.2728 


-7.182 


-0.04353 


-7.773 


-3.154 


0.686 


1.0 


7.17852 


9.1852 


-6.982 


-0.0397 


-8.393 


-3.233 


0.7099 


1.5 


6.95888 


10.863 


-7.9839 


-0.3896 


-12.29 


-4.649 


0.7491 



Table 4. Second g^-moment with respect to average, ((g^ — (g^))^) 



More detailed kinematic measurements allowing to study the dependence of the inte- 
grated rates upon correlated changes in the limits on and g°, can potentially further 
improve the sensitivity to the 'soft charm' effects by more effectively emphasizing the 
kinematics where the corresponding physics resides. A certain possibility along these 
lines has been discussed in Ref . [8] . 





V 


B 


G 


P 


D 


G 


L 





6.92486 


16.488 


-10.568 


-0.2179 


-44.39 


-12.43 


2.683 


0.5 


6.11989 


15.69 


-9.8461 


-0.2234 


-44.02 


-12.22 


2.634 


1.0 


2.60892 


10.978 


-4.9345 


-0.0397 


-44.1 


-11.35 


2.536 


1.5 


-2.03212 


7.5062 


2.9761 


-0.9282 


-62.41 


-16.01 


2.829 



Table 5. Third g^-moment with respect to average, {{q^ — {q^})^) 



Ultimately, the analysis of the inclusive semileptonic decay data should include the 
value of He in the fit, along with the other heavy quark parameters as required by the OPE. 
A simplified procedure would be using the heavy quark parameters currently extracted 
from the lepton energy and hadronic mass moments ignoring these effects, to compute 
the moments. Comparing those with the measured ones allows to infer bounds on He 
and in this way probe the nonperturbative charm contributions, or even detect them. 
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4 Discussions and conclusions 



With \Vcb\ and other heavy quark parameters hke rrib and rric having been extracted 
from B — > IvXf. with high precision, scrutinizing even small contributions comes onto 
the agenda. In this paper we have analyzed the effects of the nonperturbative QCD 
interactions of charm quarks in inclusive decays of B hadrons, specifically semileptonic 
b^ciu transitions. These effects had so far not been included in practical applications of 
the OPE; allowance for them was made in Ref. [8] in assessing the theoretical uncertainty. 

Inclusive weak decays of heavy flavor hadrons admit a local OPE, which is crucial for 
our analysis. It allows us to study nonperturbative charm effects in a model-independent 
manner. We have shown that these contributions have a well-defined meaning in the 
OPE described by expectation values of four-heavy-quark operators like {B\brccrb\B) 
(or similar operators with additional derivatives once terms of higher order in l/m^ are 
considered) . The OPE actually requires the inclusion of these four-heavy-quark operator 
terms, and they cover the impact of the soft-charm dynamics on B decays without double- 
or undercounting. For inclusive widths they scale with rrib like 

The expectation values {B\brccrb\B) can be evaluated, through an expansion in pow- 
ers of l/rric, in terms of the higher- dimensional 'usual' heavy quark operators, i.e. those 

without charm quarks. Ignoring radiative QCD effects we obtain terms scaling like 

for an overall contribution oc — Yet once hard gluons are considered, there appear 

contributions in the rate which fade out only as the first power of l/rric, as{inc) ■ In 
the picture where nonperturbative charm emerges through the higher Fock states in the 
B meson wavefunction (so-called 'Intrinsic Charm', see below) these terms describe the 
interference of the effects with and without a cc pair made possible through hard gluon 
annihilation. 

For a numerical evaluation of the nonperturbative charm effects we have elaborated 
a method for estimating the higher-dimensional b quark expectation values; it can be 
used in other applications of the heavy quark expansion as well. We also found a strong 
enhancement of the two-loop Wilson coefficients for the D = 7 operators (i.e., O 
effects) evidently related to the peculiarity of the Coulomb interaction of static quarks. 
It could potentially lead to the dominance of the O (^) contributions with individual 
ones as large as half a percent. However, we found significant cancellations between 
different contributions, see Table 2. In view of the approximations in both calculating the 
Wilson coefficients and in estimating the corresponding expectation values, we cannot take 
the resulting sub-permill numbers at face value. However, a net contribution exceeding 
0.005 Pgi looks improbable. 

As a result of the above cancellations an appreciable effect might be expected from 
D — 8 operators with the Wilson coefficients generated at one loop. The natural scale of 
these contributions can be a few permill in Fgi; explicit calculations led us to expect an 
overall effect of about 0.003 Fgi. Even considering the approximate nature of our estimates, 
we conclude that the effects associated with nonperturbative charm dynamics should not 
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downgrade the accuracy in extracting Vcb at a percent level as long as no 1/m.c expansion 
is unnecessarily built into the analysis of the semilcptonic data. 

In addition, we have suggested a more or less direct way to experimentally probe these 
effects without relying on estimates of the bbcc expectation values in a l/rric expansion. 
Such experimental tests are sensitive, of course, to the total expectation values, including 
possible exponential contributions oc e~^"*''/^'^='<*'^ which would be left out in their l/rric ex- 
pansion. The g^-distributions like the higher g^-moments are particularly constraining in 
this respect. With the present experimental capabilities it should be possible to constrain 
- or measure ~ these effects to the level relevant for the precision in Vcb down to a fraction 
of one percent. 

It has been forcefully argued in many papers over the years that hadrons may contain 
cc pairs as a higher Fock component in their wave function. This particular mechanism is 
usually referred to as 'Intrinsic Charm' (IC). More specifically it has been suggested [28] 
that such a component for B mesons could vitiate the strong CKM hierarchy in their 
decay modes expected when only 'valence' (and 'light-sea') quarks are included in the 
wavefunction. This 'Intrinsic Charm' picture can be used to visualize many of the effects 
we have discussed in this paper. However the relationship between the soft-charm effects 
discussed here and the 'Intrinsic Charm' ansatz is that of an analogy and an illustration 
rather than a genuine connection. At the technical level, the existence of a full-fledged 
OPE treatment is crucial for our analysis, and that is available for inclusive, yet not 
exclusive transitions. The intrinsic charm picture, on the other hand, taken literally 
would seem to apply to any type of process, and actually has mostly been discussed in 
the context of exclusive channels. 

Even for inclusive decays the connection is tenuous at best. The correction to Tsi{b^c) 
cannot be directly related to the admixture of cc pairs in the B meson wavefunction. The 
presence of virtual cc in the initial state does not necessarily change in a definite way the 
decay rate or other fully inclusive characteristics, much in the same way as the presence 
of the spectator quark per se does not affect the decay rate of the heavy flavor hadrons 
[1]. Charm quarks interacting with the gluon background potentially shape, to some 
extent, the structure of the B meson wavefunction and, hence, the usual nonperturbative 
parameters //^, /Iq, which affects the rates - yet the interaction of light quarks, both 
sea and valence, has conceptually the same influences. Therefore, this type of effect is not 
speciflc to 'Intrinsic Charm'. 

As discussed before, without pcrturbative corrections the leading nonperturbative 

charm expectation values are suppressed as i-e. the soft-charm correction to the 

A^ " 
width scales like . The l/m^-dependence looks consistent with a cc Fock state ad- 
mixture in the B wave function oc which leads to a term O ( ) for the cc pair 

probability. However, hard gluon corrections give rise to contributions oc ctsi^c) ^^l^i 
without an analogue in the naive probabilistic 'Intrinsic Charm' picture. 

Lastly, 'Intrinsic Charm' is meant to represent an effect due to an initial state conflg- 
uration. The OPE, on the other hand, by necessity combines initial and final state effects 
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and often does not allow a clear separation between the two. The soft extra cc pair may 
affect the decay rate showing up in the initial or in the final state, or without appearing 
at all - the OPE determines the aggregate effect of all mechanisms, including simply the 
nonperturbative corrections to the propagation of the final state quark produced in the 
decay. 

A physical interpretation of the four-quark expectation values with both b and c be- 
comes transparent if one adopts a naive factorization approximation in the usual IC 
ansatz: they would correspond to the c density at the origin, |^'c(0)P, rather than the 
overall cc-pair probability integrated over the B meson volume: 

Wcc^ J d=^x|*e-(x)|'. (66) 

In this equation we have tacitly assumed that only one pair of cc can be present, and that 
each charm antiquark in the wavefunction must be accompanied by one charm quark. 
The presence of color and spin degrees of freedom was also ignored. The four-quark 
expectation values actually select the particular projections in those spaces. 

Evidently there can be no direct relation between the admixture of a cc component in 
the wavefunction Wcc and the c density at origin |^c(0)P corresponding to the four-quark 
expectation values. Yet a general dimensional estimate 

l*c-(0)r^T^^/;cc- (67) 

involving the effective size Rq of the B meson bound state should provide some idea about 
the possible scale in relating the two quantities. They differ in dimension by mass to the 
third power, the factor which goes into the l/nif suppression of the soft-charm effects in 
the inclusive decay rates. The generic 'Intrinsic Charm' contributions do not have to be 
suppressed by powers of I /nib. According to the estimate paralleling Eqs. (43)- (46), one 
has 

_ 3 

4 / M 2 \ 2 

-TiRl ^ — ^ (0.008 GeV^)-^ (68) 
3 y 47r J 

Obviously relations of this sort cannot account for the dynamic details like the Coulomb 
enhancement we have observed in Sect. 2. 

Theory has much less to say rigorously about exclusive modes, where an 'Intrinsic 
Charm' component may naturally have a bigger impact. There arises even a fundamental 
concern inhibiting immediate answers: Does 'Intrinsic Charm' really represent a genuinely 
independent contribution to decay amplitudes when embedding parton model effects into 
real QCD, or may it in some instances involve double counting? This looks especially 
nontrivial in the context of a quantum field theory, where the highly virtual states are 
usually absorbed into various renormalization factors. Since the charm pairs in a static 
hadron have to be viewed as highly virtual, it may seem to be the case. The situation is 
additionally aggravated by the fact that the originally introduced "Penguins" as particular 
short-distance-induced local operators [29] are often used now in a rather loose and wide 
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sense. In particular, "generic Penguin" diagrams characterized only by their topology, can 
often be deformed in such a way as to seemingly exemplify the IC effects, as illustrated, 
for instance, by Figs. 1. In fact, a generic diagram for the nonperturbative correction 
to the charm propagation in the 6 — > c decays similarly includes contributions which 
can be interpreted as an IC effect with subsequent c annihilation. At the same time, it 
was shown [30] that the rescattering processes can formally 1 /mb-dominate over the bare 
annihilation-type contributions; this applies to the charm quark annihilation as well at 
sufficiently large rrih. 




Figure 1: A diagram for a 'generic Penguin' process (left) where the decay along the channel 
6 — > ccg produces a final state without hidden charm, simultaneously includes the 'Intrinsic 
Charm'-type contribution to the same final state (right) if both 'intrinsic' c and c annihilate 
with the quarks produced in the decay. 

Without going into further details here, we note that an essential stumbling block in 
clarifying these principal elements here is lack, as of now, of a consistent OPE procedure 
which could be applied to such exclusive decay processes in Minkowski space. For instance, 
a clear separation is missing of which effects are absorbed into the decay matrix elements 
and which corrections are attributed to the effective weak decay Hamiltonian. Related to 
this, there is no clear understanding of how to integrate out the high-frequency/high-mass 
modes and what type of an effective Minkowskian field theory for exclusive decays would 
result from this. A more complex answer here may, in principle, allow for more significant 
strong-interaction corrections related to non-valence charm quarks, which would be absent 
for much heavier charm quarks. 
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Appendix 

A Conventions 

Gauge fields 

The ordinary and covariant derivatives read 

da = -^, D^ = do^-ieA\ (A.l) 

with t"' — A"/2 and A" the usual hermitian SU{?>) Gell-Mann matrices 

[A", A''] = 2i/"^^A^ , Tr[A"A^] = 25'^* . (A.2) 
The field strength tensor and its dual are defined as 

Electromagnetic notations 

In terms of the field strength tensor the electromagnetic fields read 

= Foi (A.4) 

-Bj — {Fij = eijkBk) ■ (A.5) 

The square of the field strength tensor is 

F^^F^"" ^2{B^-E^) . (A.6) 

Greek indices a run from to 3, Latin indices a run from 1 to 3, and the summation 
convention is to be understood in the following way: 

a • & = a^6^ = aobo - aibi - 0262 - 03^3 = aobo - dkh ■ (A.7) 

The e-tensor 

We use the e-tensor convention from [11] which is the same as in Itzykson & Zuber as 
opposed to Bjorken & Drell, 

e'''' = -eoi23 = 1 . (A.8) 
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B Currents in the Fock-Schwinger gauge 

Here we sketch the technique used for calculating the charm currents in the external 
gluon field. We use the Fock-Schwinger gauge and work in momentum space. Our master 
equation is 

{-c'^Tc\0))a = t f ^Jh[^W^W ■ 1 , (A.9) 

where F is a matrix in Dirac space, a, h are color indices and the 'propagator' must be 
understood as an operator: 

i 1 

{k^,)ab = k^6ab - -{Ga^)ab{0)dk'- " -(^AGa,.)ab(0)<9fca(9fcA + . . . , (A.IO) 

with dk^ = We expand l/(^ — rric) in Eq. (A.9) in the gluon field strength, which 
yields a series in differential operators in the Fourier transform ('momentum') space. 
These act on the momentum space 'unity' in Eq. (A.9), which implies that the right-most 
derivative vanishes. 

To derive this representation, we start with the definition 

where ip — rric is the Dirac operator in the gluon background, and we have shown the 
spinor indices a and /? explicitly. We then write the coordinate matrix element (a;|...|a;) of 
the Dirac propagator in the momentum space: using {k\x) — e^^^ we have (the Lorentz 
indices are not shown for simplicity) 

\x) = {x\—: \x) = '^{x\p)(p\—: \k) {k\x) — 

Y^^i{j^)xi^\ ^ \k). (A.12) 



Since {p\k) = S{p—k), the latter expression upon summation over p recovers the definition 
of the operator acting on the momentum-space unity in Eq. (A.9): 

J2{p\f(k; dk)\k) = Yl f(k; dk)5(p-k) = Yl f(k; 9^) • l . (A.13) 

p,k p,k k 

The second step is to apply the Fock-Schwinger gauge {x — xo)-A{x) — 0, A{xo) — 0. 
We shall set the fixed point Xo to zero and then (for details, see Ref. ( [31]) 

A^{x) = f dXXx^G^^iXx) . (A. 14) 

t/ 

We can Taylor expand the field strength in powers of Xx and integrate over A. Due to the 
specific choice of the gauge condition the ordinary derivatives in the Taylor expansion can 
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(a) (b) (c) 

Figure 2: The one-loop graph (a) and the relevant two- loop graphs (b),(c); the external 
gauge field fines are not shown. The two-loop graph (b) is depicted in detail in Fig. 3. 

be replaced by the covariant derivatives, and the expansion of the gauge potential takes 
the following form: 

A^{x) = ^x'^G^^ - ^x'^x^D^G^^ + ... . (A.15) 

Using this explicit form for the gauge potential in the covariant derivative in the Dirac 
operator and replacing x^ by —idk„ we arrive at the representation Eqs. (A.9)-(A.10). 

C Some elements of the 2-loop calculation as{mc)/mc 

Here we outline the principal steps of the calculation that lead to the result in Eq. (21). 
The graphs are depicted in Fig 3. 

Our aim is to calculate the coefficients Cjk{ij) in Eq. (8). Using the universafity of 
the OPE we can obtain those coefficients for on-shell 6-quarks with p = {rrih, 0, 0, 0) in 
the fixed gluon background, where the matrix elements are simple functions of the gluon 
fields. 

There are six possible ways to attach a hard gluon to the one-loop graph in Fig 2a; 
only two contribute at order l/rric- The gluon has to connect the &-quark line with the c 
quark in the loop. Fig. 2b, c, otherwise the symmetry of the graph is the same as in the 
one-loop case and would bar a 1/ rric contributions. 

The sum of the two-loop graphs in Figs. 2b, c reads 

{b(p)\Jum) = (A.16) 
-g^f Tr[r2,5e(p+A;)iV>Se(g)] Aj(A;) b{p) (t''rSb{k+p)r, + r,Sb{p-k)t'r) b{p) 

Jkq ^ ' 

with Ji2 defined in Eq. (9) and the trace taken over the c-quark spinor and color indices; 
the s-quark and gluon propagators in the external field are given below, Eqs. (A.18)- 
(A.21). We use the convention Jy = J df-y and Jk = I tI^- the heavy quark limit the 
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(a) 



(b) 



(c) 






(d) 



Figure 3: The graphs corresponding to Fig. 2b which contribute to order as{m^ /itLc in the 
fixed-point gauge. The double fine shows the static h quark and the loop stands for charm 
quark. The two external gluon lines generically represent the operators in Eqs. (20), as 
it comes out naturally in this gauge. 



6-quark propagator S^, reduces to 



mi,-+oo 



ml + iO 



—^-—^±P^-m5iko). (A.17) 



Only the i/ = Q vertex survives the static limit, and the expression becomes 
{h{p)\Mh{p)) = |^^TV[r25e(^ + q^l^SM)] Koik) b {in5{ko){t', T,} - [t', T,]) b 

It can be shown that the contribution of the principal value piece vanishes upon 
integration; we have explicitly verified this in our calculation. 

As the next step we need to put the appropriate propagators in the external field and 
to project out the operators in Eq. (20). As discussed in the previous Appendix, **for 
external gluon fields we use the background field method in the Fock-Schwinger gauge 
x-A{x) = where the gauge potential is expanded directly in the gauge invariant field 
strength; c.f. Ref [31] for a pedagogical introduction. For internal gluon fields a different 
gauge may be used; we adopt the Feynman gauge for it. 

The gluonic operator parts in Eq. (20) carry dimension four and may either come from 
a single propagator, or from two separate propagators each contributing one gluonic field 
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strength of mass dimension two. It is most economic to choose the fixed point to be the 
position of the local four-quark operator. With this choice the 6-quark propagator does 
not interact with soft gluons because the Wilson line is unity in this gauge; other soft 
contributions to the b quark propagator vanish in the heavy quark limit. There remain 
six graphs shown in Fig. 3 out of ten possible combinatorial possibilities for Fig. 2b, plus 
similar diagrams corresponding to Fig. 2c. They will be examined below. 



Propagators in the external field 

The massless fermion propagator in the external field was given in Ref. [31], Eqs. (2.30)- 
(2.31);^ the gauge potential propagator in the external field was also presented in [31] 
for the case where one of the endpoints is the fixed point, Eq. (2.35). For the fermion 
propagator we need to include mass. The diagrams in Fig. 3 also include gluon Green 
function between two general points. Both generalizations are not difficult. 

We arrange the propagators in increasing dimension of the external gluon fields: 

(x, y)= -t {TAl{x)AUy)) = A^^^\x, y) + Ag]'^''(a;, y) + Ag]" ''(x, y) + ^^^l^^x, y) 

S{x,y) = -i{T^{x)^{y))A^S^'^\x,y) + S^''\x,y)+S^^\x,y) + S^^\x,y), (A.18) 

The propagators are obtained using the background field Lagrangian. With the fixed 
point a;o = the used terms for the fermion propagator are 



S^'\x,y) 


= S{x-y) = / 


S^'\x,y) 


= -^G./3(0)| 


S^'\x,y) 









{idp + x)a 
{-idq + y)a 



e-'P<^^S{p)^f,S{q) , 
{y-idq)a{y-idq)fj{y-idq)j i Jp 



idp x)q,( idpj^q x)^ 
-idq+k + y)a{-idk + y)^ 
id„ + x)Jid„ + x)g(id„ + x 



and for gluons 

-1 



A2^''(a;, y) = g,, 5^' e^^^^"^) — , ( A.20) 

Ag]«''(x,y) = -2GJ(0) / 

+ ^^..GS(-^4 + |/)a/^e-H^)pl, 



g ip{x-^)]_ 



^The DDG term in S^^") was omitted in [31], Eq. (2.31), since the authors aimed at calculating only 
the gluon condensate contribution (G^). 
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jp p 



+ 



For the fermion propagators we gave different possible representations which can be chosen 
pursuant to calculational convenience, see below. The momenta q and k must be set equal 
to p after taking derivatives. For the gluon fields the matrix notation = f"'^'^G^^^ is 
assumed. In fact, one can do without a coordinate representation if one of the endpoints 
of the propagator coincides with the fixed point. This is the case for the c quark lines, 
but not for the gluon propagator in Fig. 3. S^^^ or A*^^), respectively, do not contribute, 
and we do not give 5"*^^^ and A^"^^ for brevity. 

After putting in the propagators and evaluating the traces the needed integral takes 
the generic form 

Those integrals are in fact easily calculated in arbitrary dimensions d. Rewriting the 
propagators as {q^—rn^Y = (— l)"(g^+(m^— gp))'' suggests integrating over d'^^q and d'^^k 
first. We combine the denominators by introducing a Feynman-Schwinger parameter 
which is integrated over at the end yielding the Euler Beta-function. The variable ko is 
fixed by the (5-function and the remaining integration over qo is elementary. 
For the first step the following master integral is needed: 



J(a, b, c;d) — J 



d'^k d'^q 1 



{27rY {27rY {k^Yiq^+mYiik + qY+m^Y 

2' ' 2' 



T{^ — a)T{a+h+c—d) , . ... os^„jja^-i 



{^^Y r(6)r(c)r(f) 

and, therefore the scalar integral in Eq. (A.21) is 

d'^fc d'^g dgo 



/oo(a, c;d+l) = y 



— — ' , , , —B{a^h-%,a^c-%){m^Y^'^^. (A.23) 

[^^Y^\ r(6)r(c)r(f) ^ 2' ^ ^ ^ 

In view of proliferating combinatorial possibilities we implemented the evaluation of the 
tensor integrals algorithmically starting from the scalar master integral Eq. (A.23). 
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The c-quark loop potentially has UV divergences and the gluon loop potentially has 
IR divergences. Since the one-loop contribution to this order in l/nic vanishes, the di- 
vergences must disappear in the sum over all the graphs. The only diagrams involving 
IR or UV divergent integrals are those from Figs. 3 d,e,f which are those where a gluon is 
emitted from the hard gluon line. Explicit calculations yield that, in fact, the divergences 
cancel in each individual graph. 

Furthermore, a Furry-type analysis of the symmetries of the graphs can be done. Since 
the Gell-Mann matrices have no definite transformation properties under transposition, 
one has to resort to Hermitian conjugation. This distinguishes the QCD analysis from 
the QED case. The outcome is that the two graphs in Fig. 3d,e are equal; the two graphs 
in Fig. 3a,b are equal in the vector channel, but not in the axial channel. This is verified 
explicitly in the calculations. Another nontrivial consistency check would be to introduce 
a generic Fock-Schwinger fixed point and to observe its disappearance from the final result. 
This would require, however more involved calculations and the evaluation of four more 
graphs. 
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